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Abstract 



O I We apply expansion methods to obtain an approximate expression in terms 

[ of elementary functions for the space and time dependence of wave packets 

^~»' in a dispersive medium. The specific application to pulses in a cold plasma is 

considered in detail, and the explicit analytic formula that results is provided. 
When certain general initial conditions are satisfied, these expressions describe 
^SJ ' the packet evolution quite well. We conclude by employing the method to 

o ■ 

' exhibit aspects of dispersive pulse propagation in a cold plasma, and suggest 

' how predicted and experimental effects may be compared to improve the 

Q ' theoretical description of a medium's dispersive properties. 
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I. INTRODUCTION 



A systematic study of the propagation of electromagnetic pulses in a dispersive medium 
originated with Sommerfeld [|T],3 and Brillouin The subject has subsequently been 

studied by means of a number of interesting methods and novel interpretations; see for 
example references and references therein. Some of these methods are extensions to 
the classic Sommerfeld-Brillouin theory. 

The purpose of this paper is firstly to provide an analytic expression for the behavior 
of dispersive wave packets by use of a straightforward method [see Eq. (|])], along with a 
description of the initial conditions that must be satisfied in order to render this expression 
valid. Secondly, we furnish the explicit formula that results when the method is applied 
to the case of Gaussian pulses in a cold plasma [see Eq. ([T6|)]. We conclude by using 
the method to discuss and exhibit features of the dispersive propagation and interaction of 
multiple pulses in such a plasma. 



II. GENERAL METHOD 

The method we employ is based on Taylor expansion techniques and Gaussian wave 
packet expansion techniques that allow required integrals to be computed exactly. More 
specifically, we expand an initial wave packet as a superposition of Gaussian packets; then 
power-expand the dispersion relation to second order around the dominant wavenumber 
of each Gaussian. This allows for an exact solution of the Fourier integral describing the 
evolution of the packet. 

We take as our starting point the standard linear Fourier integral [0,|rn[ describing the 
propagation of a pulse E{x,t) in a medium of dispersion uj{k): 

poo 

E{x,t)= A{k)e'^^''-^^''^'Uk, (1) 



where A{k) must be chosen to satisfy the initial boundary conditions for E{x,t) For 
simplicity, we assume that -^E{x, 0) = 0, so that A{k) is just the Fourier transform of the 
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initial pulse E{x, 0) = Eq{x) P,[T0 



1 

A{k) = — EQ{x)e-'^'' dx. (2) 

We then take the real part of Eq. (|1]), assume A{k) is real, and denote Re£'(x, t) by E{x, t) 
for simplicity 

/oo 
y4(A;) cos (fcx - uj{k)t) dk. (3) 

The next step is to find a fit to Eq{x) in the form of a superposition of sinusoidal wave 
packets with a Gaussian envelope 



n 



Eo{x) = J2fie-'^^^ cos(M, (4) 
1=1 

where the constants a^, bi, fi are real, and Oj > 0. Such a fit can normally be found for a 
localized initial packet Eq{x). The Fourier transform of Eq. (|) can then be found exactly 

as 



1 /-oo 1 " f- 

A{k) = - Eo{x)e-^'^ dx = -^Y: 



(5) 



This amplitude is real, so the assumption below Eq. is valid. The dominant wave 
numbers in Eq. (^) are ±6j, i = l,n, so we expand the dispersion relation uj{k) to second 
order around each point ±6j, i = l,n, separately: 

u{k, ±bi) = So{±bi) + ks,i±k) + kh2{±k). (6) 

Equation (j^) constitutes 2n expansions of uj{k), each expansion being centered around one of 
the 2n points ±6j, i = l,n. The coefficients Sj(±6j) will in general depend on the parameters 
describing the dispersion relation uj{k) of the dispersive medium. 

We can now obtain a good approximation to Eq. by using Eqs. (|^) and (|^) to expand 
the integrand in Eq. (0) as a sum of 2n terms, where each term is centered on the major 
wave numbers ±bi, i = l,n of the initial packet: 
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E{x,t) 



1 " f- 



4v^ •=! ^ 



-tso{bi) + k{x - tsi{bi)) - tk'^S2{bi 



dk 



(7) 



Note that Eq. (|^ arises from a generalized expansion approach to Eq. since it incor- 
porates an arbitrary number of (i2j^erent expansion points ±6i. The integrand in Eq. is 
expanded to a sum of 2n Gaussian terms by substituting Eq. (||) into Eq. (||), and the only 
modification of a term occurs via a Taylor expansion (|]) of the factor uj{k) around the point 
where the term is most influential. The advantage of the expansion (|^) is that the resulting 
integrals can now be found analytically as 

1 " f { 
E{x,t) = -J2^\F^t,a,,b,) 



X exp \F{t, Gi, 6i)(4aj 



[x - tsi{bi)] tbiS2{bi) - ^ [a; - tsi(bi)f - t^blsKpi 



X cos Oj, bi 



{Aai)-^{[x - tsi{h)] bi - t [bls2ih) + soibi)] } 



+ -[x - tsi{k)f ts2{bi) - thl{k)so{bi 



--tan [AtaiS2{bi) 



(8) 



where F(t,ai,bi) = [(Aai)'"^ + 1'^ sl{bi)] Equation (H) is useful in the sense that it is 
an explicit formula containing only elementary functions, and in the sense that it is quite 
general. All one needs to do to apply it to a specific situation is to compute the power series 
of the dispersion relation around the main wave numbers of the initial wave packet. The 
accuracy of the expression (P) will of course increase as the number of expansion points are 
increased. The formula provides for expansions of uj{k) up to second order, but one could 
also of course limit the order to unity by taking S2(±&i) = in Eq. (^. 

We now discuss the initial conditions that must hold for Eq. (|^) to be valid. From Eq. 
(^), we see that the amplitudes of the two terms of the integrand in the contribution 

1 /. 
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COS {kx — Lu{k)t) dk 



(9) 



to E{x,t) in Eq. (j^) are maximal at the points k = bi and k = —hi. The amphtudes decay 
exponentially as k moves away from the points ±6j, and are of their maximal value 
when is at a distance 2^/5^ from either point. At a distance S^/oi, the amplitudes are 
g-6.25 ^ ]^ g ^ io~3 Qf their maximal value, which is so small that beyond this distance one 
may consider the contribution (|^) to be negligible. The /c-intervals 

I{hi) = [bi - by/al, bi + S^/o^, I{-bi) = [-bi - by/al, -bi + (10) 

may therefore be taken as the only two areas of the fc-line where the contribution @ to 
is appreciable. 

In order for Eq. (H) to be applicable, one must require that a contribution (for a 
specific i) is appreciable only within the two areas of the fc-line where the two second order 
expansions uj{k, ±6j) of uj{k) (for the above i) in Eq. (^ are valid. The remainders \uj{k) — 
uj{k, bi)\ and \uj{k) — uj{k, —bi)\ will in general be small over a certain interval containing bi 
and another interval containing {—bi), respectively. We denote the largest two such intervals 
over which \uj{k) —uj{k, and \uj{k) —uj{k, —bi) \ are small by I'{bi) and I'{—bi) respectively. 
These intervals depend on uj{k) and ±6j [The remainder terms depend on the third derivative 
of uj{k), as is well known from Taylor expansion theory]. The requirement for Eq. (^ to be 
valid can accordingly be stated as 

/(6,) C I\h), C I'{-bi), (z = l,n). (11) 

Therefore, when the intervals /'(&«) and I'{—bi) have been determined from uj{k) and 
±6j, the domain of validity of Eq. (^ can be stated as an upper limit on the pulse shape 
parameters in Eq. (^, as noted from Eq. ([T0|). Note that the n criteria in Eq. (pUD 
are independent, since the resultant packet E{x,t) in Eq. (|^) is a linear combination of n 
independent contributions of the form (|^), as seen from Eq. (|^). 
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III. COLD PLASMA APPLICATION 



We now consider the specific case of a pulse propagating in a cold plasma. For simplicity, 
we assume that the initial pulse can be approximated well by a single normalized Gaussian 
packet centered at x = 0, 

Eo(x) = e-'^^'cos(te). (12) 

In this case, Eq. must be applied with the parameters 

n = 1, = a, bi = b, fi = 1. (13) 

For the cold plasma, we use the dispersion relation [|r^ 

Uj{k) = ^c^k^ + cu2, (14) 

where c is the speed of light in vacuum, and Up is the plasma frequency. 

We then expand Eq. ([T^) around b and —b, to obtain the coefficients in Eq. (P) as 



so{-b) = so{b), si{-b) = -si{b), S2{~b) = S2{b), (15) 



where B = {b'^c^ + ujp)^^- These coefficients are general, since they are functions of the 
initial pulse parameter b. 

Substitution of Eqs. (|^) and (|15D in Eq. (||) yields 



E{x, t) = ^a-^Ti exp{(x - tb^c^B^tbc^ujl{Sa)~^B^ - {x - tb^c^B^f{16a)'^T 
-t^b'c^uj^^{16a)-'B'} 

X cos{T[— (x - tb^c''Bl)a-% + -{x - tb^^Blftc^ulB'^ 
16 8 ^ 

- Ua-^B'i{2b\^ujl + ul) - h'u;yBl{3bW^ + 2u;'^)] 

- I tan-\2tac^ujlBl)} + {b ^ -b), (16) 
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where T = (^a~^ + jt'^c^ujpB^)~^ and B = {b'^c^ + ujp)~^. Equation (|T6D is a general 
equation, describing the dispersive behavior of pulses of the form ([T^ that propagate in a 
medium with dispersion of the form (p!^). 

The intervals /(±6j) and the validity criteria for Eq. (0) are found by substitution of 
Eq. (p!3D in Eqs. (0) and (p^. Since the dispersion relation uj{k) in Eq. ( |T^ ) is symmetric 
about A; = 0, the validity criteria in Eq. ([iTl) simplify. A second order expansion uj{k, b) of 



the square root uj{k) = ^c^k"^ + uj^ about a positive wave number b approximates uoik) quite 
well when > 0, but fails when A; < 0. One must therefore require the interval I{b) in Eq. 
(|T0|) to lie to the right of the origin of the fc-line. According to Eq. (0), the left endpoint 
of I{b) is to the right of the origin when b > 5^/a. 1(b) then lies within [0, b + 5^/a\, which 
is within the interval I'{b) of validity of the expansion uj{k,b) [given by Eqs. (^) and (|15|)] 



of u{k) [given by Eq. (0)]. A similar result holds for the expansion u{k, —b) of ^c^k"^ + oj^ 
about a negative wave number —b. Hence, we have the useful initial condition that 

\ > 5 (17) 



for Eq. to be valid. In other words, the approximation (|TB]) is accurate when the 

frequency b of the initial pulse Eq{x) in Eq. (|T^ is high, or when Eq{x) is broad (i.e. when 
the pulse shape parameter a is small). However, we see from Eq. ( [TtI) that short pulses are 
also well approximated by Eq. (|16D when they are of high frequency, as well as low-frequency 
pulses of long duration. 

In the specific illustrations and discussions below, the pulse parameters a and b were of 
course always chosen so that Eq. (^) was satisfied, with many choices such that ^ >> 5. 
As a double check, we also compared values for E{x,t) given by Eq. (|1^) with the values 
for E{x,t) obtained by numerically integrating Eq. [using Eqs. (|), (|13), ( [1^ ) and (^^] 
for a systematic variation of values for x, t, a, b, c and Up. For values of a and b obeying 
Eq. (p!7D, the agreement was excellent, as expected. The discrepancy between the exact and 
approximate values for E{x,t) was typically found to be less than 10~^. 



Equation ([Iq ) represents two pulses traveling in opposite directions, as we expect flO 



from the initial condition 4i?(a;, 0) = 0, stated above Eq. @. For simplicity, in the 



remaining part of this paper we focus only on the right-traveling pulse term in Eq. ([TBD 
(which is the first term in Eq. (|16D when b is positive). 

We now outline some of the features of the dispersive behavior of a single pulse in a 
cold plasma, which can be observed from plots of Eq. (|1^). It is readily observed that 
the envelope of a short initial pulse broadens more rapidly as time elapses compared to the 
broadening of the envelope of a long initial pulse, the initial pulse frequencies being equal. 
This is of course well known, and can be explained analytically . 

It is also seen that in general, the number of oscillations within the pulse envelope in- 
creases with time, and more interestingly, the effect of dispersion on the pulse is to shift high- 
frequency components of the pulse toward the spatial front of the pulse and low-frequency 
components of the pulse toward the spatial rear of the pulse (see the pulse El in Fig. H). 
This effect becomes more prominent as the initial pulse is shortened. 

Additionally, the spatial packet velocity of a pulse increases as the pulse frequency is 
increased (with spatial packet velocity we mean Xmax/t, where Xmax is the distance traveled 
by the spatial maximum of the pulse envelope during a given time period t. Temporal 
pulse velocity is similarly defined as x/tmax, where tmax is the time that elapses before the 
temporal maximum of the pulse envelope appears at a given distance x [^]). However, as 
expected, in no cases is the center of a packet beyond the point x = ct aX time t if it is at 
the origin x = at time t = 0. Here c is the velocity of light in vacuum, appearing in Eq. 
(|T6|) from the dispersion relation (p!^ . 

Moreover, the envelope shape is better maintained for a high frequency initial pulse 
than for a low frequency initial pulse, the initial pulse durations being equal. Therefore, 
increasing the frequency of a short pulse will improve the preservation of its shape as it 
travels through a dispersive medium, as well as increase its packet speed. The propagation 
of short, high-frequency pulses is a topic of strong current interest [|12[, [1^ . 



We next examine the the dispersive evolution of more than one (right-traveling) pulse 
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through a cold plasma. One can readily do this by repetitive use of Eq. (|Tq). For example, 
the propagation of two initial pulses -^01(2;) = e~'*i^^ cos(6ix) and -^02(2;) = e""^^^^ cos(62a^) 
separated by an initial time delay d is described by 

E{x,t) = El{x,t) + E2{x,t), (18) 

where El{x,t) and E2{x,t) are obtained by substituting (t — t + c/, a — > ai, 6 — > bi) and 
(a ^ a2, b ^ 62) respectively, into Eq. ([TBI). The first pulse El is here d time units ahead of 
the second pulse E2 when E2 passes the origin x = at the time t = {A description of the 
general dispersive behavior of m initial pulses Eqj of the form (^ separated by initial time 
delays dj can analogously be obtained from Eq. (H) by the set of substitutions t t + dj, 
Qi aij, bi — > bij for i = l,n and j = l,m ). 

Figures |l] through ^ depict some features of dispersive double pulse propagation, as 
described by Eqs. (|16D and ([T8|). These figures also illustrate some of the properties of 
single pulse propagation mentioned above (note that the packets in the top plots of Figs. |1| 
and 1^ at the initial time t = are identical to those given by Eq. (|T^), except that their 
amplitudes have half the magnitude of those in Eq. (|T2|). This is because half of the energy 



of the initial packets is carried by the left-traveling parts of E{x, t), which are omitted from 
all plots). 

Figure shows how two initially distinct packets of the same frequency and duration 
gradually overlap and interact with each other due to the dispersive stretching of the two 
packets as time elapses. Figure ^ shows the results of the interaction after a very long time. 
One can clearly see the emergence of constructive and destructive interference effects among 
the different spatial regions of the coalescing pulses. 

In Fig. ^, a high frequency, short pulse is dispatched after the dispatch of a low frequency, 
long pulse. We see that the high frequency pulse eventually overtakes the low frequency 
pulse, in accord with the single-pulse feature that the packet velocity is greater for a high- 
frequency initial pulse than a low-frequency initial pulse. The middle and bottom plots show 
how the high frequency pulse interferes with the various spatial parts of the low frequency 
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pulse. The bottom plot depicts an instant of constructive interference between the two 
pulses, where the total amplitude is larger than the separate amplitudes of either pulse. For 
longer times than those shown in the figure, the high frequency pulse eventually passes and 
leaves the other pulse. 

As a final example of a usage of Eq. (|TB|), one may employ this Eq. to determine 
parameters a^, bj, and initial time delays dj, {j = i,m) of m pulses such that they will all 
have the same width and completely overlap at a specified time t while propagating in a 
plasma of dispersion characterized by c, Up and Eq. (0) . Figure ^ shows a two-pulse version 
of this, with t = 1000. We see that strong destructive and constructive interference effects 
occur in this case. A comparison of such predicted effects with those found experimentally 
in a specific medium could be used to determine the medium's dispersion relation u!{k). 

IV. SUMMARY AND CONCLUSIONS 

We have provided an analytic expression describing the propagation of dispersive wave 
packets [Eq. (|^)] provided the packets satisfy the applicability criteria (pT|). The expression 
is obtained by a clear-cut method, and can be used to study properties of the propagation 
process (for example pulse velocity and multiple pulse interference effects due to dispersion). 

In the specific case of packets propagating in a cold plasma, we used Eq. (||) with two 
terms, and provided the explicit, analytic expression that results in this case [Eqs. (|l^) and 
(0)]- For initial pulses satisfying the condition (|17D , Eq. (|16D gives an accurate description 
of single or multiple pulse propagation in a dispersive cold plasma. 

As a final note, it is conceivable that the presence or non-presence of effects predicted 
by Eqs. (P) and ([16|) could be used to verify, improve the equations for, or determine a 
substance's dispersive properties. 
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FIGURES 

FIG. 1. Temporal evolution of two pulses El (right) and E2 (left) in a cold plasma, obtained 
from Eqs. ( [T^ ) and ([l8|). The initial time delay between the two packets is d = 1.5. Both pulses are 
identical, with initial- value parameters ai = a2 = 10 and 6i = 62 = 20^. The plasma is character- 
ized by c = 1 and ujp = 1. The spatial distribution of the total field E{x,t) = El{x,t) + E2{x,t) 
is shown at times t = (top), t = 500 (middle) and t = 1000 (bottom). 

FIG. 2. Same as in Fig. |, but for t = 5000 (top) and t = 100000 (bottom). 

FIG. 3. Temporal evolution of two pulses El (right) and E2 (left) in a plasma, where El is 
characterized by ai = 5, 61 = 20^, E2 by 02 = 100, 62 = 60^, and the plasma by c = 1, = 1 
and Eq. (0). -El is initially d = 1.1 time units ahead of E2. The spatial distribution of the total 
field E{x,t) = El{x,t) + E2{x,t) is shown at times t = (top), t = 500 (middle) and t = 1000 
(bottom). 

FIG. 4. Spatial dependence of a packet El{x, t) (top) and a packet E2{x, t) (middle) after they 
have both traversed a cold plasma (characterized by c = 1, Up = 1 and Eq. ([T^)) for a period 
of t = 1000 time units. Initially, El was d = 1.1 time units ahead of E2. El is characterized 
by ai = 10, bi = 20^ and E2 by 02 = 2, 62 = 40^^. The bottom plot shows the total field 
E{x, 1000) = El{x, 1000) + E2ix, 1000). 
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FIGURE 1 
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FIGURE 2 
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FIGURE 3 
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FIGURE 4 
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